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Abstract 

We investigate the potential of a flavor pure high gamma electron capture electron 
neutrino beam directed towards a large Water Cherenkov detector with 500 kt fiducial 
mass. The energy of the neutrinos is reconstructed by the position measurement within 
the detector and superb energy resolution capabilities could be achieved. We estimate 
the requirements for such a scenario to be competitive to a neutrino/ anti- neutrino 
running at a neutrino factory with less accurate energy resolution. Although the 
requirements turn out to be extreme, in principle such a scenario could achieve as 
good abilities to resolve correlations and degeneracies in the search for sin 2 2#i3 and 
5cp as a standard neutrino factory experiment. 
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1 Introduction 



All observations on neutrinos coming from the sun [1-10], the atmosphere [11-18], and re- 
actors [19, 20] are well understood in the picture of neutrino oscillations [21] in the three 
generation framework of lepton mixing. Two of the mixing angles, sin 2 2# 12 and sin 2 2# 2 3 
have been measured as well as the two mass square differences (Am 2 -,] and Am 2 x have been 
determined. Furthermore, the parameters which are mainly relevant in the atmospheric 
neutrino oscillations, i.e. sin 2 2# 2 3 and jAm^l have been confirmed by the terrestrial exper- 
iments K2K [22-24] and MINOS [25,26]. 

However, the remaining two mixing parameters, the third mixing angle sin 2 2# 13 and the CP 
violating phase <5cp have not been determined yet. Currently, there only exists an upper 
bound for sin 2 2# 13 [27,28] and there is no information on the value of 5 C p- Also, the sign 
of the mass squared difference Am 2 x is currently unknown, i.e. it is unclear if neutrinos 
exist in normal or inverted hierarchy. So, the aim of future oscillation experiments is to 
measure these two parameters, to improve the precision to the leading solar and atmospheric 
parameters, and determine the neutrino mass hierarchy. In order to complete the picture of 
neutrino oscillation parameters, several types of new experiments have been proposed and 
are studied extensively. This includes reactor experiments [29-34] that are only sensitive to 
sin 2 2#i3, and experiments where information on both, sin 2 2^13 and 5qp can be obtained, 
like superbeam experiments [35-40], neutrino factories [41-49], and beta-beams [50-63]. 

Recently, another idea has been proposed, which makes use of a neutrino beam with neutri- 
nos coming from electron capture processes [64,65]. The electron neutrinos that are emitted 
from such electron capture processes would have a definite energy Q in the rest frame of 
the mother nuclei. Therefore by accelerating the mother nuclei to a Lorentz factor 7 the 
neutrino energy E u can be completely controlled, since the energy of the neutrinos that are 
boosted exactly towards the direction of the detector is E u = 2 / ~fQ. So, the 7 factor and 
the baseline length L have to be chosen respectively to the Q value of the electron capture 
process, the location of the oscillation maximum, and the minimal energy observable at the 
detector, e.g. above the Cherenkov threshold of muons at a Water Cherenkov detector. For 
example, if Q is relatively large (0(1 MeV)), 7 can be chosen to be of the order (9(100). In 
this case the neutrino beam can be viewed as exactly monoenergetic in the detector [63-70]. 
On the contrary, if Q is relatively small ((9(100 keV)) the 7 must be chosen quite high 
(0(1000)), but the necessary choice of the baseline leads to the effect that the neutrinos 
now have a wider energy range at the detector. While the maximal energy of E u = 2 / ~fQ 
is reached by the neutrinos in the beam axis, the energy of the neutrinos becomes smaller 
off the axis and the minimal observable energy of the neutrinos depends on the detector 
size and the baseline. In this scenario, the neutrino energy can be reconstructed from the 
vertex position measurement relatively to the beam axis within the detector and in principle 
a superb energy resolution can be achieved [64,66]. This however requires in addition to 
the resolution of the position measurement within the detector, that the beam divergence 
of the stored mother nuclei can be accurately controlled. This scenario seems interesting 
since only with one acceleration factor 7 a wide range of neutrino energy can be covered 
simultaneously with a very accurate neutrino energy determination. 

In this work we investigate the potential of such scenarios with a flavor pure electron neutrino 
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beam coming from beta capture at high 7 lead towards a Water Cherenkov detector with 
a fiducial mass of 500 kt. We will refer to these scenarios as monobeam scenarios in the 
following. We estimate the requirements for such a scenario to be able to resolve correlations 
and degeneracies in the search for the remaining oscillation parameters sin 2 2# 13 and Sep 
within the measurement in only one polarity, i. e. neutrino running, but with superb energy 
resolution abilities and to be competitive to a standard neutrino factory scenario with 
neutrino and anti-neutrino running, but less accurate energy reconstruction. Unfortunately, 
the ability to also gain information on the sign of Am^ at the discussed monobeam scenarios 
is limited due to the missing anti-neutrino running, so it will be omitted throughout this 
work. 

This work is organized as follows: In Section 2 we sketch the basic principles of the high 
gamma electron capture monobeam experiments and summarize all underlying assumptions. 
Furthermore, we define the reference setups that are investigated throughout the rest of 
the work. Next, in Section 3 we address the issue of requirements to resolve parameter 
correlations and degeneracies in the sensitivity to sin 2 2#i3 at the reference scenarios defined 
in Section 2 and then address the sensitivity to 5qp in Section 4. Here, also all parameter 
correlations and degeneracies are taken into account. We summarize and conclude the 
main results in Section 5. Finally, the details of the operation of a monobeam experiment 
including the energy reconstruction by the position measurement, i.e. a derivation of the 
neutrino energy E U (R) as a function of the radius from the beam axis, and the details of 
the event rate calculation is presented in the Appendix. 

2 Experiment configurations and simulation techniques 

In the scope of this work we consider a flavor pure neutrino beam that is produced within 
the electron capture process of ^Sn isotopes: 



In the rest frame of the process the produced neutrinos are monochromatic with an energy 
of Q = 267 ke\0 at a lifetime of 4.11 h_|. The isotopes are assumed to be accelerated in 
a decay ring, where they coincide with electrons accelerated at the same 7 factor and a 
boosted neutrino beam is produced towards the detector. At the distance of the baseline L 
the neutrinos hit the detector at a radial distance R from the beam axis and their energy 

1 We only consider electron capture from the K shell here. A more detailed analysis should also include 
electron capture from higher shells. However, the results should not be affected significantly. On one hand a 
position measurement of a neutrino would allow different true energy values and a new discrete uncertainty 
arises, but on the other hand the ratio is known and the higher the shell, the smaller the contribution. 
Furthermore, the sets of neutrinos from electron captures from other shells cannot be interpreted as back- 
ground since also their energy is accurately known, besides a discrete uncertainty, and they also oscillate 
and hence contribute to the fit. 

2 This is the main limiting factor for obtaining an adequate number of electron capture processes per 
year, i.e. to collect enough statistics. However, in [71-73] the possibility to enhance the electron capture 
rate has been discussed. 
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in the laboratory frame (rest frame of the detector) can be expressed as 
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The derivation of this formula can be found in the Appendix, whereas the approximation is 
taken from [64] and can be obtained in the limit of large 7 with (3 ~ 1 — ^2 and (R/L) <C 1. 
At the beam center the neutrino energy is maximal at a value of E u = and decreases 
for larger distances from the beam center. Since the neutrino energy is a function of the 
distance from the beam center, a position measurement within the detector allows a precise 
reconstruction of the neutrino energy. We assume a Water Cherenkov detector with a 
fiducial mass of 500 kt. The large detector mass allows to collect enough statistics that 
is needed to gain from the superb energy resolution and can have large geometrical size 
in order to have a enough broad energy window, since the minimal measurable energy 
depends on the maximal distance from the beam center. We assume the geometry of the 
detector to be as shown in Figure [TJ The radius of the detector is set to -R max = 100 m 
so that the depth is still approximately 64 m and a reconstruction of the Cherenkov rings 
and electron/muon identification remains possible. The position measurement should be 
optimized for this kind of experimental setup and reach at least a resolution of AR =30 cm, 
which has been the estimated vertex resolution at Super-K for fully-contained single ring 
events [14]. Furthermore, the vertex resolution for muon events, i.e. the monobeam signal 
events, is slightly better than for electron events and can even reach a resolution of 25 cm 
in the energy window of interest. It should be mentioned that the very good position 
measurement resolution can only be transfered into an excellent energy resolution if the 
systematical uncertainty in the beam spread can be reasonably controlled. This means 
that the divergence of the stored isotopes perpendicular to the beam line must satisfy the 
condition p x /p z < AR/L before the decay. Otherwise, the superb energy resolution that is 
assumed in this work could not be achieved although the position measurement is accurate. 
This translates for baselines that are discussed in the following into the requirement of a 
beam divergence p x /p z < 1 yurad for the mother nuclei in the storage ring and seems hardly 
feasible. However, it should be noted that beam divergences of 0(1 fir ad) are already 
discussed, for instance for the proton beam of the LHC for the operation of the TOTEM 
experiment [74]. 

For having neutrino energies beyond the Cherenkov threshold and allow for electron/muon 
discrimination, we only discuss monobeam setups with neutrino energies above 400 MeV. 
The signal efficiencies and background rejection factors follow the analysis from [54] (mainly 
the low gamma beta beam therein). Above 400 MeV up to 700 MeV there was found a 
signal efficiency of approximately 0.55 for the appearance measurement of muon neutrinos, 
which we take to be the signal efficiency of the discussed monobeam scenarios. Although the 
signal efficiency in [54] decreases slightly for higher neutrino energies, we assume the signal 
efficiency to stay stable up to the highest energies discussed in this work (Ej, < 1.4GeV), 
since the monobeam setup does not rely on the quasi-elastic events only, because the energy 
reconstruction can be performed by the position measurement within the detector. The 
rejection factors for background coming from neutral current events with single pion pro- 
duction at energies above 400 MeV are found to be below 10~ 3 in [54], whereas we assume 
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Figure 1: Scheme of the detector setting. The fiducial volume is indicated by the dashed lines. The 
neutrino beam hits the detector at the edge of the fiducial volume to allow for distance measurements from 
the beam axis. In case of i? max = 100 m the depth of the fiducial volume part would be approximately 64 m 
for a 500 kt fiducial detector mass. 

this background rejection to be at a level of 10~ 4 . This is a quite optimistic assumption, 
and it is not clear, that this could be achieved. However, note that the mismatch of ordi- 
nary energy reconstruction and energy reconstruction by position measurement due to the 
carried away missing energy of the neutrino in neutral current reactions could give a further 
rejection ability of such background events. We assume a systematical uncertainty of 2.5% 
for the signal events and 5% for the background events, as also assumed for the typical 
beta beam scenarios, i.e. in [61]. The uncertainty of the signal events has probably to be 
called optimistic, but since we will find that the main effect will come from correlations and 
degeneracies [75-77], a value of 5% would not have much impact to the results of this work. 

As can be understood from Eq. (T5|), the energy window of the analysis is, due to the technique 
of energy reconstruction, limited by the size of the detector to the interval 



,., 2 l Q /n2 <E V < 2 7 Q, (3) 
1 + (7 J R max /L) 2 

so that the energy window is completely fixed after the baseline L and the acceleration 
factor of the ions 7 is chosen. So finding an optimal Setup is more complicated as it is 
for example in the case of beta beams, since choosing a perfect pair of L and 7 to exactly 
measure at the first oscillation maximum can suffer from an energy window that is to small 
to allow resolving correlations and degeneracies. However, adjusting the baseline to smaller 
baselines in order to have a lower minimal energy also shifts the oscillation maximum to 
lower energies, while going to higher values of 7 not only shifts the maximal energy but also 
the minimal energy to higher values. So, the whole energy window moves away from the 
oscillation maximum although it is broadened. Therefore, in the next sections we discuss 
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the potential and performance of the following different reference scenarios of monobeam 
setups: 

• Setup I: The Water Cherenkov detector with a fiducial mass of 500 kt is located at 
a baseline of L=600km, the mother nuclei 5Q°Sn are accelerated with 7 = 2500 and 10 
years of data taking are assumed at the number of 10 18 electron capture decays per 
year. 

• Setup II: The Water Cherenkov detector with a fiducial mass of 500 kt is located at 
a baseline of L=250km, the mother nuclei 5Q°Sn are accelerated with 7 = 2000 and 10 
years of data taking are assumed at the number of 10 18 electron capture decays per 
year 

• Setup III: The Water Cherenkov detector with a fiducial mass of 500 kt is located 
at a baseline of L=600km, the mother nuclei 5Q°Sn are accelerated with 7 = 900 and 
7 = 2500 consecutively, and 5 years of data taking are assumed in each of the two 
phases so that as for Setup I and II the total running time is 10 years. The number 
of 10 18 electron capture decays per year is assumed for both phases. 

Setup I is located at the first oscillation maximum, but the energy window is not very 
broad compared to the width oft the oscillation maximum peak, therefore we also discuss 
the second scenario, Setup II, with a broader energy window which on the other hand is 
located slightly off the first oscillation maximum at higher neutrino energies due to the 
smaller baseline. Then again, because of the smaller baseline higher event rates can be 
obtained at Setup II. With Setup III we discuss the potential in resolving the correlations 
and degeneracies with a monobeam experiment by a combination of data from the first 
oscillation maximum and also the second oscillation maximum. This combination should 
be a powerful tool to resolve the degeneracies and the importance of the second oscillation 
maximum has been discussed in [78] . Since the first oscillation maximum phase at Setup III 
is comparable to Setup I, the gain from the additional measurement at the second oscillation 
maximum can directly be read off the comparison of Setup I and Setup III. The exact 
width of the corresponding energy windows of the setups and their location respectively 
to the oscillation maxima are shown in Figure [2j Note, that Setup III makes use of the 
combination of different 7 which was also the strategy for the "high Q-low 7" electron 
capture beam experiment scenarios as discussed in [65,67-70]. However, for these scenarios 
this strategy was required to obtain spectral information at the first oscillation maximum, 
while Setup III provides spectral information at the higher 7 = 2500 and information from 
the second oscillation maximum is included with the second arrangement of 7 = 900. This 
can be seen in Figure [2j 

There, the appearance probability P(u e — > is plotted for sin 2 2^13 = 0.01 and three 
choices of #cp (the other oscillation parameters are chosen as in Eq. (j3])). The yellow (grey) 
bands indicate the energy window of the analysis for Setup I and III in the left-hand side 
and Setup II in the right-hand side. It can be seen that the energy window for the choice 
of L=600km and 7 = 900 is essentially only a very narrow band while for the higher values 
of 7 indeed a broader energy window can be covered over the whole radius of the detector. 
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Figure 2: The appearance probability P(z^ e — ► z^) as a function of the neutrino energy E„ at a baseline 
of L=250km (left-hand side) and L=600km (right-hand side). The oscillation parameter values are the ones 
from Eq. ([!]), sin 2 29 = 0.01 and the values for the phase 5cp are chosen as labeled in the plot legend. The 
vertical yellow (grey) bands indicate the energy window of the analysis that is given for a detector radius of 
Rmax = 100m, the respective baseline L and the chosen 7 factors, i.e. 7 = 9OO/7 = 2500 for the scenarios 
at L=600km (Setup I and Setup III) and 7 = 2000 for the scenario at L=250km (Setup II). 

However, the energy window of Setup I is too narrow to cover the first oscillation maximum 
for the different choices of 5cp- For 5qp = the peak of the first oscillation maximum lies 
inside the energy window of the analysis but for the maximally CP violating values for Sep 
the peak moves outside the energy window. The energy window of Setup II lies above the 
first oscillation maximum independent of 5qp but we will show in the next sections that 
Setup I will suffer more from correlations and degeneracies than Setup II since the latter 
benefits from a higher event rate due to the smaller baseline and the larger energy window, 
where the superb energy resolution can evolve. 

Note, that the number of electron capture decays per year taken for the reference scenarios 
is of the order of the "high Q" electron capture scenarios discussed in [65,67-70] and also 
the order of beta decays per year discussed for the beta beam scenarios [50-62] . However, 
because of the long lifetime of the 5o°Sn electron capture of 4.11 h in the rest frame the 
feasibility to achieve a number of 10 18 per year seems out of range if it cannot be enhanced. 
This enhancement of the electron capture processes could be achieved by a laser irradiation 
as discussed in [71-73]. Furthermore, as is also the case for high gamma beta beams [54,61] 
the high gamma values require a very large accelerator complex of the size of the LHC. 

For reasons of comparison and to put the performance of the monobeam setups into per- 
spective we will compare the results to a standard neutrino factory setup with a 50 kt MID 
detector at a baseline of L=3000km and a parent energy of the stored muons = 50GeV. 
This neutrino factory setup is similar to the standard scenario for a neutrino experiment [45], 
commonly known as NuFact-II, with 1.06- 10 21 useful muon decays per year (corresponding 
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to 5.3- 10 20 useful muon decays per year and polarity for a simultaneous operation with both 
polarities). The details of the neutrino factory experiment description follow the description 
of the NuFact-II scenario in [49] . We assume a runtime of five years in each polarity, so that 
the total running time is 10 years as for the discussed monobeam setups. Furthermore we 
will also consider an optimized neutrino factory scenario at the end of Section 4, where com- 
pared to the standard neutrino factory scenario a second detector similar to the standard 
detector at L=3000km is installed at the magic baseline L=7500kmjf] 

The analysis throughout this work is performed with the GLoBES software [80, 81] and 
the incorporated Poisson x 2 -analysis. Details can be found in the GLoBES manual [82]. 
Since the monobeam only measures ^-appearance and could additionally only observe v e - 
disappearance, the leading atmospheric parameters sin 2 2023 and jAm^l cannot be deter- 
mined as would be the case at a neutrino factory with a measurement in the z/^-disappearance 
channel. Thus, correlations with the leading atmospheric parameters would spoil the po- 
tential of the monobeam experiment alone, as also would be the case for a beta beam 
for the same reasons. Therfore, we adopt the same technique as in [61] and add the v^- 
disappearance information from a simulation of the superbeam experiment T2K. The cor- 
responding appearance information is excluded, so that information on sin 2 2#i3 and #cp is 
solely collected by the monobeam experiment (see [61] for details). The errors on the solar 
parameters are taken to be 5% on each, Am^ and #i 2 . 

As input or so-called true values within the simulations, we use, unless stated otherwise the 
following parameter values, close to the current best fit values (see Refs. [83-86]): 

Am 2 x = 2.5 • 10~ 3 eV 2 sin 2 2 9 23 = 1.0 , 

Am 2 21 = 8.2 • 10~ 5 eV 2 sin 2 26 12 = 0.83. (4) 

Note, that the octant-degeneracy [87] does not affect our results, since we choose sin 2 2# 23 
to be maximal and thus the octant-degenerate solution appears at the same point in the 
parameter space as the original solution in the parameter space. So, if it is stated that 
effects of degeneracies are taken into account in the analysis, only the intrinsic sin 2 29^- 
#cp-degeneracy [88] and the sign(Am| 1 )-degeneracy [76] are regarded out of the complete 
set of the so-called eight-fold degeneracy [77]. 

3 Sensitivity to sin 2 2#i 3 

The sensitivity to sin 2 2# 13 is calculated under the hypothesis of true sin 2 26 , 13 = 0. The 
sensitivity limit at a certain confidence level is then the maximal fit value of sin 2 26 Vi that 
still fits the simulated data at the chosen confidence level, i.e. it would be the lower bound 
to sin 2 2^13 that the experiment could achieve in case of vanishing true sin 2 2^13. It is 
well known, that the main problem is to resolve the correlations with the other oscillation 
parameters and the so-called eight-fold degeneracy. In Figure [3] the sensitivity to sin 2 2# 13 
is shown at the 3cr confidence level as a function of the number of decays per year for 

3 The optimized scenario furthermore uses an optimized disappearance channel with the MINOS energy 
threshold [38] while the muon CID with the implied CID cut threshold is only used for the golden appearance 
channel. See [79] for details. 
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Figure 3: The sensitivity to sin 2 26 1% at the 3a confidence level for the monobeam scenarios 
L=600km/7 = 2500, L=250km/7 = 2000, and L=600km/ 7 = 900 + 7 = 2500 as a function of the num- 
ber of decaying ions per year including statistics, systematics, correlations, and degeneracies. The lowest 
curve represents the pure statistical sensitivity limit to sin 2 2#i3 and the colored bands indicate the effect of 
switching on systematics (blue/dark grey), correlations (green/middle grey), and degneracies (yellow/bright 
grey) so that the final sensitivity limit is given by the upper curve. 

the monobeam scenarios at L=600km/7 = 2500, L=250km/7 = 2000, and L=600km/7 = 
900 + 7 = 2500. The vertical lines indicate the reference setups at a number of 10 18 ion 
decays per year. In each plot the lowest curve represents the pure statistical limit to 6 
and the colored bands show how the sensitivity degrades if also systematics (blue/dark 
grey band), correlations (green/middle grey band), and degeneracies (yellow/bright grey 
band) are taken into account. The final achievable sensitivity limit to sin 2 2#i3 is given 
by the upper curve. Obviously the statistical and systematical sensitivity limit to sin 2 2# 13 
at all three scenarios in Figure [3] can reach to very small values of sin 2 2# 13 due to the 
very large statistics in the Water Cherenkov detector. However, the monobeam scenario at 
L=600km/7 = 2500 can resolve the correlations not until an exposure of 10 17 decays per 
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year. The point where the degeneracies can be resolved is reached not until approximately 
10 20 decays per year, which of course is beyond any feasibility. So despite the improvement 
of the statistical limit with higher exposures the final sensitivity limit to sin 2 29 13 stays 
relatively stable a approximately sin 2 26 13 m independent of the number of decays per 
year. The monobeam scenario at a baseline of L=250km and 7 = 2000 suffers from the 
same problem. First, the sensitivity limit does only slightly improve and almost stays stable. 
Beyond exposures of 10 18 decays per year this scenario starts to resolve the degeneracies and 
the sensitivity limit to sin 2 2#i 3 improves significantly. From Figure [3] it becomes clear, that 
the technique of a high gamma monobeam with its superb energy resolution in a narrow 
energy window is not able to resolve the correlations and degeneracies in a measurement 
at just one 7. The scenario at a baseline of L=600km allows to measure in the second 
oscillation maximum since for L=600km this maximum is located above the Cherenkov 
threshold and events can be collected. The lower plot of Figure [3] shows the sensitivity limit 
to sin 2 2# 13 for such a scenario, where 5 years data taking at 7 = 900 and 5 years data taking 
at 7 = 2500 is combined. Now, the correlations and degeneracies can be already resolved 
for lower exposures. We checked that it is not necessary to split up the two data taking 
phases into an equal period of five years each. The ability to resolve the correlations and 
degeneracies still remains if only 2 years data taking at 7 = 900 are combined with 8 years 
at 7 = 2500 and the final sensitivity would be even slightly better since then more statistics 
could be collected at the first oscillation maximum. 

For reasons of comparison, the sensitivity to sin 2 28 13 at Setup I, Setup II, and Setup III 
are again shown in the left-hand side of Figure H] and confronted with the sensitivity limit 
obtainable at the standard neutrino factory scenario. The neutrino factory also suffers from 
the correlations and degeneracies. But as can be seen in the right-hand side of Figure H] the 
difference is that the neutrino factory can almost resolve the degenerate solution. There, 
the projected A% 2 is shown as a function of the fit value of sin 2 2# 13 for the degenerate 
solution with the wrong sign, i.e. inverted hierarchy while the positive Arn^j was taken as 
input true value. The degenerate solution appears for the neutrino factory scenario at a A% 2 
only slightly below the 3cr, while the degenerate solution for Setup I appears at A% 2 = 
and thus fits as good as sin 2 2# 13 = 0. On the other hand, with Setup III there does 
not appear a second local minimum in the projected A% 2 so the combination of first and 
second oscillation maximum data gives a strong tool to resolve the degeneracy. However, 
resolving the degeneracies remains the main problem if one want to reach to very small 
values of sin 2 2#i 3 and one could also think of a combination of a monobeam setups with 
the anti-neutrino running of a standard beta beam scenario. It should be noted that the 
performance of a neutrino factory could be improved by additional data from the silver 
channel — > u T [89,90], a second detector at the magic baseline [77,91,92] or a lower 
threshold (see [79]). 

4 Sensitivity to CP violation 

Due to the continuous intrinsic sin 2 26 l i 3 -5cp-degeneracy a total rates analysis of appearance 
data of neutrinos only would give continuous bands as allowed regions in the sin 2 2#i 3 - 
<5cp plane. If combined with a second band from appearance data of anti-neutrinos only 
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Figure 4: Left-hand side: Comparison of the sensitivity to sin 2 26*i3 at the 3a confidence level at the 
monobeam scenarios Setup I, Setup II, Setup II, and a neutrino factory at a baseline of L=3000km including 
systematics, correlations and degeneracies. The left edge of the bars indicates a pure statistical sensitivity 
limit. The right edges of the bars indicate the sensitivity limit after switching on systematics (blue/dark 
grey), correlations (green/middle grey), and correlations (yellow/bright grey), so that the rightmost edge 
gives the final sensitivity limit to sin 2 26*13. Right-hand side: The projected A% 2 as a function of the fit 
value of sin 2 2t9 13 fitted under the assumption of inverted hierarchy while the true values are given with 
sin 2 2#i3 and normal hierarchy. The rightmost intersections of the curves with the grey horizontal 3a line 
give the right edges of the bars in the plot on the left-hand side. 



two intersections, the true and the degenerate allowed region remain. Adding the spectral 
information obtained with conventional energy resolution, the degenerate solution can be 
resolved in most cases. This is the planned procedure at superbeam experiments, neutrino 
factories as well as beta beam experiments to resolve the sin 2 2# 13 -(!>cp-degeneracy. However, 
at a monobeam experiment only neutrino appearance is observable and the question arises, 
if and under which circumstances the superb energy resolution abilities of a monobeam 
could in principle compete in resolving the sin 2 26 , 13 -c J >cp-degeneracy. Since we found in 
the last section that the ability in resolving the degeneracies does not appear until a large 
number of decays per year, we will fix this value to 10 18 decays per year in all the following 
considerations and only discuss the fixed scenarios Setup I, Setup II, and Setup II. In 
Figure [5] the allowed regions in the sin 2 2f9 13 -<5cp-plane at the 3a confidence level are shown 
for different choices of input true values. This figure is for illustrative purposes only and 
no correlations with the other oscillation parameters is considered, i.e. they are kept fixed 
to the values of Eq. (j4]). The left column is for Setup I (L=600km/7 = 2500), the middle 
column is for the Setup II (L=250km/7 = 2000), and the right column shows the allowed 
regions obtained for the standard neutrino factory setup for reasons of comparison. The 
bands, indicated by the solid grey lines, represent the corresponding allowed regions at 
the 3<T confidence level if only total rates are considered. As expected, the total rates 
allowed regions for the monobeam scenarios are bands that do not restrict <5cp at all whereas 
for the neutrino factory already also the parameter space of Sep is restricted due to the 
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Figure 5: The allowed regions in the sin 2 2#i3-<$cp-plane for the true values indicated by the black dots 
at 3<7 for Setup I (left column), Setup II (middle column), and a standard neutrino factory (right column) 
for reasons of comparison. Only the correlation between sin 2 29±3 and 5qp are taken into account and all 
other parameters are fixed to values of Eq. ((4]) . The plots also contain the allowed regions at 3<r for total 
rates only (grey solid lines). 
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Figure 6: Sensitivity to any CP violation at 1 (yellow/bright grey), 2, 3, 4, and 5a (red/dark grey) 
after 10 years of data taking as a function of the true values of sin 2 29i3 an d Sqp- The sensitivities are 
shown for the monobeam scenarios Setup I (upper left-hand side plot), Setup I (upper right-hand side plot), 
Setup III (lower left-hand side plot) and a standard neutrino factory (lower right-hand side plot) for reasons 
of comparison. For a pair of true values within the shaded regions the CP conserving fit values <5cp = 
and (5cp = 7r can be excluded at the respective confidence level. 

information from neutrino and anti-neutrino data. If spectral information is included to the 
analysis, the neutrino factory allowed regions are not influenced significantly and only the 
small degenerate solutions can be excluded, but for the monobeam scenarios because of the 
superb energy resolution wide parts of the bands can be excluded and only smaller allowed 
regions remain that are comparable in size to the allowed regions from the neutrino factory 
scenario. However, in some cases of choices of true values still degenerate solutions remain. 
As mentioned before, we have ignored correlations with the other oscillation parameters and 
also the sign(Am| 1 )-degeneracy here. In all of the further considerations, we will focus on 
the sensitivity to CP violation if also these correlations and all degeneracies are taken into 
account. 
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The sensitivity to any CP violation is shown in Figure [6] for Setup I (upper left-hand side 
plot), Setup II (upper right-hand side plot), Setup II (lower left-hand side plot), and the 
neutrino factory scenario (lower right-hand side plot) at the 1, 2, 3, 4, and 5 a confidence 
level from bright grey/yellow (lex) to red/dark grey (5a). Sensitivity to any CP violation 
is given for a pair of true values sin 2 2#i3-<5cp if the CP conserving values <5cp = and 
Sep = 7r do not fit the simulated reference data if all correlations and degeneracies are taken 
into account. It is known, that the standard neutrino factory suffers from the sign(Amg 1 )- 
degeneracy in some areas of the parameter space (sin 2 29 13 w 10~ 2 ' 5 and 5qp ~ —n/2), 
because of the so-called "7r-transit" , i.e. the degenerate solution fitted with wrong sign of 
Anijj contains the CP conserving value for 8qp = vr (see [49] for details). As can be seen from 
Figure El Setup I suffers strongly from correlations and degeneracies at larger true values of 
sin 2 26*13 whereas Setup II performs better. Within the interval <5cp £ [— tt, 0] Setup II does 
not suffer from any correlations and degeneracies anymore and gives better results than the 
neutrino factory in the same interval. In the interval Sqp £ [0, 7r] Setup II and the neutrino 
factory perform in a comparable manner, only for larger true values of sin 2 2#i3 > 10 -2 the 
neutrino factory looses sensitivity to CP violation for values of <5cp near the CP conserving 
values. This effect is due to the uncertainty of the matter density along the baseline which 
strongly affects the performance of a neutrino factory at large values of sin 2 2# 13 because 
of the very long baseline. The best sensitivity to any CP violation is found for Setup III. 
Here, the combination of data from the first and second oscillation maximum can resolve the 
degeneracies that appear at the baseline of L=600km for Setup I. Additionally the sensitivity 
to CP violation of Setup III reaches to significant smaller values of sin 2 26i% at the maximally 
CP violating values <5cp = ±7r/2. We checked that, as also was the case for sensitivity to 
sin 2 2^ 13 , a combination of 2 years at 7 = 900 and 8 years at 7 = 2500 would also already 
allow to give this performance. The results from Figure are finally summarized in Figure [3 
The fraction of <5cp parameter space where sensitivity to any CP violation is given at the 3a 
confidence level is shown as a function of true sin 2 2#i3 for the considered scenarios Setup I, 
Setup II, Setup III and a neutrino factory. The performance of the standard neutrino factory 
is indicated by the black solid line. However, we also show the performance of an optimized 
neutrino factory scenario, where in addition to the standard golden channel measurements 
at the baseline L ~ 3000km a second 50kt Magnetized detector is installed approximately 
at the magic baseline L = 7500km. As can be seen in Figure [7J the performance of the 
neutrino factory setup is significantly improved. Note, that a CP fraction of 1 can never 
be achieved, since values near the CP conserving values can never be distinguished due to 
finite statistics. 

5 Summary and conclusions 

We have analyzed the potential of high gamma neutrino beams from electron capture decays 
of 50° Sn isotopes directed towards a large Water Cherenkov detector with a fiducial mass 
of 500 kt. The resulting neutrino beam would be completely flavor pure and only consist 
of electron neutrinos. The achievable resolution in the energy reconstruction in such a 
scenario can be significantly more precise than from the usual energy reconstruction in Water 
Cherenkov detectors, since it is performed by the position measurement within the detector. 
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CP violation 




Figure 7: Comparison of the fraction of (true) <5cp for which CP violation can be established at the 3<r 
confidence level as a function of the true value of sin 2 2#i3 at the monobeam scenarios Setup I, Setup II, 
Setup III. The solid black line is for a standard neutrino factory while the dashed line is for an optimized 
neutrino factory with a second detector at the magic baseline. 

The aim of this work was to estimate the potential and requirements of such scenarios to 
resolve the correlations and degeneracies in the sensitivity to sin 2 2#i3 and the sensitivity 
to any CP violation, only with their power in energy resolution abilities. This has been 
compared to the performance at a neutrino factory, where the combination of neutrino- and 
anti-neutrino running is used to resolve correlations and degeneracies. We have compared 
three monobeam setups, two of them with a different energy window at different locations 
respective to the first oscillation maximum. Setup I at a baseline of L=600km and 7 = 2500 
has been chosen such, that the energy window of the analysis is directly located at the 
first oscillation maximum, but due to this choice the energy window is not broad enough 
to cover the whole oscillation maximum. Setup II at a baseline of L=250km and 7 = 2000 
on the other hand has a broader energy window which is located at higher energies as the 
oscillation maximum. In comparison to Setup I this setup gains from the broader energy 
window and the larger statistics due to the smaller baseline. The baseline of L=600km also 
allows to take data at the second oscillation maximum, which is at this baseline already 
located at energies above the Cherenkov threshold of muons. Therefore Setup III combines 
a measurement at the first oscillation maximum (7 = 2500 as in Setup I) and the second 
oscillation maximum (7 = 900), 5 years data taking each. For the exposure of the setups it 
has been assumed to have a running time of 10 years at a number of 10 18 decays per year. 
This number is hard to obtain because of the relative long lifetime of the 5Q°Sn isotopes and 
an enhancement of the electron capture rate has to be achieved. However, concerning the 
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sensitivity to sin 2 2#i3 we found that this number is required to evolve capabilities to start 
resolving the correlations but still the pure superb energy resolution and the high statistics 
alone cannot compete with the sensitivity to sin 2 29 13 at a standard neutrino factory with a 
50 kt MID detector at a baseline of L=3000km and a parent muon energy of E M = 50GeV 
because the degeneracies cannot completely be resolved. On the other hand the neutrino 
factory also suffers strongly from degeneracies and additional data from the silver channel, 
the magic baseline or lower energies (maybe with a hybrid detector) would be required. 
Setup III on the other hand with the combination of data from first and second oscillation 
maximum performs well in resolving the correlations and degeneracies. It gives a better 
sensitivity sin 2 2#i3 < 2.5 • 10~ 4 at the 3a confidence level. When it comes to the sensitivity 
to any CP violation the performance of the monobeam setups is more impressive than the 
performance concerning the sensitivity to sin 2 2# 13 . Setup I still suffers significantly from the 
degeneracies while Setup II reaches sensitivity in a quite large part of the parameter space 
and no negative effects from degeneracies could be observed. Finally, Setup III showed very 
good abilities to establish CP violation in a very large part of the parameter space and all 
degeneracies coming from the measurement in Setup I can be resolved due to the data from 
the second oscillation maximum although no information from a anti-neutrino running is 
included. However, one has to note that the requirements to achieve such a performance, i.e. 
the very large acceleration factors of the isotopes, the high number of isotope decays per 
year, and the very low beam divergence of the stored isotopes of (9(1 /irad) are extreme. 
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A Relativistic transformations 



A.l Neutrino energy 

In the following considerations, the primed quantities are defined as the quantities in the 
laboratory frame, i.e. the rest frame of the detector and the quantities without a prime 
represent those in the rest frame of the electron capture decays in which the neutrinos are 
produced. 



Rest Frame 



Lab Frame 





In the rest frame of the decays the neutrinos are produced at an energy Q and with an 
uniform angular distribution of momentum. Since for the considered mother nuclei ^Sn 
the endpoint energy is Q = 267 keV (m„ <C Q), the neutrino mass can be neglected: 

Q 2 =p 2 + m 2 ^p 2 . (5) 

So, for a neutrino that is emitted in the direction r the four-momentum in the rest frame 
of the decay is given by 

( Q \ ( Q \ 



v » = 



V J 



Q sin 9 cos</> 
Q sin 9 sin (f> 
Q cos 9 



(6) 



Since the problem is 0-symmetric, we can choose = for the sake of simplicity and the 
four-momentum of the neutrino in the rest frame can be written as 



P 



(7) 



/ Q " 

Q sin 9 


\ Q cos 9 ) 

The boost is in the z-direction, and after the transformation the energy of the neutrino in 
the lab frame becomes 



E' = p'° = 7 Q(l +f3 cos9). 



(8) 



A. 2 Transformation of angles 



Now we want to derive the energy of a neutrino that hits the detector at a baseline L and 
at the distance R from the beam center, i.e. at an angle 

L 1 



cos#' 



VL 2 + R 2 yj 1 + (L/Kj 



(9) 
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The expression for the neutrino energy has to be found as a function of the angle cos 6', or 
respectively the radius R. 

From p^' it is quite straight forward to find the transformation of cos 9: 

7Q(/3 + cos6>) (3 + cos 9 



cos 9' = — = ... 

v / [tW + cos 9)f + (Q sin 9) 2 l + f3cos9 

and in the other direction the transformation is given by 

. -/3 + costf' 

cos 9 = — . 

1 -/3cos0' 

The transformation of is trivial = 0' and therefore we find that 

d£l d cos 9 



dVL' dcos9' 



with 



and the corresponding 



dcos9 



dcos9 



d cos 9' 
dcos9 



-= [ 7 2 (l-/?cos£') 2 ] 



21 -1 



= [ 7 2 (l+/3cos0) 2 ] 



21-1 



(10) 



(11) 



(12) 



(13) 



(14) 



Now, the exact formula for the neutrino energy in the lab frame as a function of the lab 
frame quantities is found to be: 



£7'(cos0') 



Q 



7 1 - 13 cos 9' 



(15) 



and 



E\R) 



Q 

7 



1 - 



P 



A/1 + WLf 



(16) 



B Calculation of event rates 



The initial neutrino beam consists only of electron neutrinos. In the detector the muon 
neutrinos from the appearance channel will be detected. The neutrino energy is maximal 
at the beam center (R = 0) with E max = 2 / ~fQ and decreases to the outer regions of the 
detector. We introduce an equidistant binning in R 2 to have more balanced event numbers 
in the different bins, than would be the case for equidistant binning in R. In the simulations, 
we use k = 100 bins, so that the largest bin appears in the beam center with approximately 
10 m radius and the smallest bin is found at the outer edge of the detector with a width 
of approximately 50 cm, whereas the position measurement resolution is assumed to be at 
least 30 cm, which is the vertex resolution estimated for fully-contained single ring events 
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at Super-K [14] in the energy window of interest. In this work we have not introduced an 
additional smearing between the bins in the outer regions of the detector. However, if the 
vertex resolution cannot be optimized this binning turns out to be too narrow at the outer 
bins in the detector and smearing would have to be introduced to these bins or the width 
of the bins in the analysis would have to be customized. We checked, that going to an 
equidistant binning in R 2 with only 50 bins, i.e. bin sizes from 14 m to 1 m, or going to an 
equidistant binning in R with 100 bins hardly changes the main results of this work. 

For the usage within the GLoBES software, the radial binning is translated to binning in 
energy, where the bins are not equidistant anymore. 

If -R max is divided in k bins the edges of the bins are 

= RL, -(i- i)A^ 2 (17) 

with 

d2 

AR 2 = (18) 
k 

Here always R 2 > R 2 + i holds, so that the corresponding energy bins are in the right ordering 
for GLoBES : 

E\R 2 ) < E\R 2 +l ). (19) 
Furthermore, within GLoBES for the calculations the mean value of each energy bin is taken: 

E, = g(fi?) + g( ^>. (20) 
2 

Then, the appearance event number in one energy-bin is given by 

I on 

Ni-ax P(L, Ei)^ x ^-(£?) x a{E>) x N QUC>i , (21) 

where q is the signal efficiency in the corresponding bin, P(L, Ei) Ve ^ v is the appearance 
oscillation probability, is the angular neutrino flux, a(E' i ) is the charged current 

cross section per nucleon, and Af nuCi i is the number of nucleons within the geometrical size 
of the i-th bin: 



Here is the fraction of all number of nucleons that have to be considered in the i-th energy 
bin. 

Since the neutrino flux in the rest frame of the decays is uniformly distributed, it can be 
written as 

d^l -^decays {01\ 
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where iVdecays is just the number of decays, i.e. the number of produced neutrinos. The 
neutrino flux can now be found with Eqs. ([13]) and (fT5l) : 



dn_ dn dtt _ iV decays ■■ 2 fl^fl/^T 1 _ decays ( Ej V 

dtl'i dQdQ' Arc 171 P cosu *>\ 4n \Q J ' 

Also, it is straight forward to show by using Eqs. (|T2l) . (fl3|) . ffT5"]) , and (|23|) that 

(in (ifi cZO' A^ecays 



(24) 



dn = — ——dE' = t^UdE', (25) 
dQdQ'dE' 2/? 7 Q ' v ; 



i.e. is constant. 



References 



[1 
[2 

[3] 
[4] 
[5] 

[6] 

[7 
[8 
[9 
[10 

[11 
[12 

[13 

[14 
[15 



B. T. Cleveland et al, Astrophys. J. 496, 505 (1998). 

J. N. Abdurashitov et al. (SAGE), J. Exp. Theor. Phys. 95, 181 (2002), 
|astro-ph/0204245 , 



W. Hampel et al. (GALLEX), Phys. Lett. B447, 127 (1999). 



M. Altmann et al. (GNO), Phys. Lett. B490, 16 (2000), [hepjjex/0006034 



S. Fukuda et al. (Super-Kamiokande), Phys. Rev. Lett. 86, 5656 (2001), 



hep-ex/010 3033, 



S. Fukuda et al. (Super-Kamiokande), Phys. Rev. Lett. 86, 5651 (2001), 
|hep-ex/0103032| 

S. Fukuda et al. (Super-Kamiokande), Phys. Lett. B539, 179 (2002), |hep z ex/0205075l 



Q. R. Ahmad et al. (SNO), Phys. Rev. Lett. 89, 011302 (2002), |iuicl r e x/0204009 



Q. R. Ahmad et al. (SNO), Phys. Rev. Lett. 89, 011301 (2002), |^cT^7p204008 



S. N. Ahmed et al. (SNO), Phys. Rev. Lett. 92, 181301 (2004), [nucl :: e x/0309004. 

Y. Fukuda et al. (Super-Kamiokande), Phys. Rev. Lett. 81, 1562 (1998), 
|hep-ex/9807003[ 

Y. Fukuda et al. (Super-Kamiokande), Phys. Rev. Lett. 82, 2644 (1999), 
hep-ex/9812014[ 



S. Fukuda et al. (Super-Kamiokande Collab.), Phys. Rev. Lett. 85, 3999 (2000), 



hep-ex/0009001 



Y. Ashie et al. (Super-Kamiokande), Phys. Rev. D71, 112005 (2005), hep-ex/0501064 



M. Ambrosio et al. (MACRO Collab.), Phys. Lett. B434, 451 (1998), [5^6x79807005 

19 



[16] 
[17] 
[18] 
[19] 
[20] 
[21] 
[22] 
[23] 
[24] 
[25] 
[26] 
[27] 
[28] 
[29] 

[30] 

[31] 
[32] 
[33] 
[34] 

[35] 
[36] 
[37] 
[38] 
[39] 
[40] 



F. Ronga, Nucl. Phys. Proc. Suppl. 100, 113 (2001). 



M. C. Sanchez et al. (Soudan 2), Phys. Rev. D68, 113004 (2003), hep-ex/0307069 
M. Ambrosio et al. (MACRO), Eur. Phys. J. C36, 323 (2004). 



K. Eguchi et al. (KamLAND), Phys. Rev. Lett. 90, 021802 (2003), hep-ex/0212021 



T. Araki et al. (KamLAND), Phys. Rev. Lett. 94, 081801 (2005), Eep ex/0406035 



Z. Maki, M. Nakagawa, and S. Sakata, Prog. Theor. Phys. 28, 870 (1962). 



M. H. Ahn et al. (K2K), Phys. Rev. Lett. 90, 041801 (2003), [hep 1 ex/0212 007 
Y. Oyama (2005), [hip 1 ex/05 12041] 



M. H. Ahn et al. (K2K), Phys. Rev. D74, 072003 (2006), hep- ex/0606032 



N. Tagg et al. (MINOS) (2006), |hep2ex70605058 , 

D. G. Michael et al. (MINOS), Phys. Rev. Lett. 97, 191801 (2006), |h1^x70607088 , 



M. Apollonio et al. (CHOOZ), Phys. Lett. B466, 415 (1999), hep-ex/9907037 



M. Apollonio et al. (CHOOZ), Eur. Phys. J. C27, 331 (2003), hep-ex/0301017 



H. Minakata, H. Sugiyama, O. Yasuda, K. Inoue, and F. Suekane, Phys. Rev. D68, 



033017 (2003), hep-ph/ 021111lj 



P. Huber, M. Lindner, T. Schwetz, and W. Winter, Nucl. Phys. B665, 487 (2003), 
hep-ph/0303232l 



K. Anderson et al. (2004), hep-ex/0402041 



F. Ardellier et al. (2004), hep-ex/0405032 



F. Ardellier et al. (Double Chooz) (2006), [hep~-ex/0606025| 

P. Huber, J. Kopp, M. Lindner, M. Rolinec, and W. Winter, JHEP 05, 072 (2006), 
|hep-ph/0601266[ 



H. Minakata and H. Nunokawa, Phys. Lett. B495, 369 (2000), |hep~-ph/0004114j 
J. Sato, Nucl. Instrum. Meth. A472, 434 (2001), [hep~-ph/000 8056. 



B. Richter (2000), hep-ph/0008222 



E. Abies et al. (MINOS) FERMILAB-PROPOSAL-P-875. 

Y. Itow et al, Nucl. Phys. Proc. Suppl. Ill, 146 (2001), |hep£ex70106019} 

I. Ambats et al. (NOvA) FERMILAB-PROPOSAL-0929. 



20 



[41] S. Geer, Phys. Rev. D57, 6989 (1998), hep-ph/9712290 



[42] A. De Rujula, M. B. Gavela, and R Hernandez, Nucl. Phys. B547, 21 (1999), 



hep-ph/9811390. 



[43] K. Dick, M. Freund, M. Lindner, and A. Romanino, Nucl. Phys. B562, 29 (1999), 
|hep-ph/9903308 



[44] V. D. Barger, S. Geer, and K. Whisnant, Phys. Rev. D61, 053004 (2000), 
|hep-ph/9906487, 



[45] A. Cervera et a/., Nucl. Phys. B579, 17 (2000), hep-ph/0002108 



[46] C. Albright et al, Nucl. Phys. B547, 21 (2000), hep- ex/0008064 



[47] A. Blondel et al, Nucl. Instrum. Meth. A451, 102 (2000). 



[48] M. Apollonio et al. (2002), hep-ph/0 210192 



[49] P. Huber, M. Lindner, and W. Winter, Nucl. Phys. B645, 3 (2002), |hep-p h/0204352 
[50] P. Zucchelli, Phys. Lett. B532, 166 (2002). 



[51] M. Mezzetto, J. Phys. G29, 1771 (2003), hep-ex/0302007 



[52] J. Bouchez, M. Lindroos, and M. Mezzetto, AIP Conf. Proc. 721, 37 (2004), 
hep-ex/0310059l 



[53] M. Mezzetto, Nucl. Phys. Proc. Suppl. 143, 309 (2005), hep-ex/0410083 



[54] J. Burguet-Castell, D. Casper, J. J. Gomez-Cadenas, P. Hernandez, and F. Sanchez, 
Nucl. Phys. B695, 217 (2004), ^ep r ph/0312068| 

[55] J. Burguet-Castell, D. Casper, E. Couce, J. J. Gomez-Cadenas, and P. Hernandez, 



Nucl. Phys. B725, 306 (2005), hep-ph/0503021 



[56] F. Terranova, A. Marotta, P. Migliozzi, and M. Spinetti, Eur. Phys. J. C38, 69 (2004), 



hep-ph/0405081. 



[57] A. Donini, E. Fernandez- Martinez, P. Migliozzi, S. Rigolin, and L. Scotto Lavina, Nucl. 



Phys. B710, 402 (2005), hep-ph/0406132 



[58] A. Donini, E. Fernandez-Martinez, and S. Rigolin, Phys. Lett. B621, 276 (2005), 



hep-ph/0411402 



[59] A. Donini and E. Fernandez-Martinez, Phys. Lett. B641, 432 (2006), [hep"-ph/060326l| 



[60] A. Donini et al. (2006), [hip z ph/0604229 



[61] P. Huber, M. Lindner, M. Rolinec, and W. Winter, Phys. Rev. D73, 053002 (2006), 
|hep-ph/0506237[ 



21 



[62] J. E. Campagne, M. Maltoni, M. Mezzetto, and T. Schwetz (2006), hep-ph/0603172 



[63] C. Volpe, J. Phys. G34, Rl (2007), 15^^/0605033 



[64] J. Sato, Phys. Rev. Lett. 95, 131804 (2005), hep-ph/0503144 



[65] J. Bernabeu, J. Burguet-Castell, C. Espinoza, and M. Lindroos, JHEP 12, 014 (2005), 
|hep-ph/0505054 



[66] J. Sato, Nucl. Phys. Proc. Suppl. 155, 180 (2006). 

[67] M. Lindroos, J. Bernabeu, J. Burguet-Castell, and C. Espinoza, PoS HEP2005, 365 
(2006). 

[68] J. Bernabeu, J. Burguet-Castell, C. Espinoza, and M. Lindroos (2005), 
|hep-ph/0512299[ 

[69] J. Bernabeu, J. Burguet-Castell, and C. Espinoza, Pos HEP2005, 182 (2006), 
|hep-ph/0512297[ 

[70] J. Bernabeu, J. Burguet-Castell, C. Espinoza, and M. Lindroos, Nucl. Phys. Proc. 
Suppl. 155, 222 (2006), [hlp^ph/0510278| 

[71] M. Ikeda, I. Nakano, M. Sakuda, R. Tanaka, and M. Yoshimura (2005), 



hep-ph/0506062. 



[72] M. Yoshimura (2005), [hep 1 ph/0507248| 



[73] T. Nomura, J. Sato, and T. Shimomura (2006), hep-ph/0605031 



[74] M. Deile (TOTEM) (2004), hep-ex/0410084 



[75] M. Koike, T. Ota, and J. Sato, Phys. Rev. D65, 053015 (2002), hep-ph/0011387 



[76] H. Minakata and H. Nunokawa, JHEP 10, 001 (2001), hep-ph/0108085 



[77] V. Barger, D. Marfatia, and K. Whisnant, Phys. Rev. D65, 073023 (2002), 



hep-ph/0112119 



[78] J. Arafune, M. Koike, and J. Sato, Phys. Rev. D56, 3093 (1997), hep-ph/9703351 



[79] P. Hub er, M. Lin dner, M. Rolinec, and W. Winter, Phys. Rev. D74, 073003 (2006), 
|hep-ph/0606119[ 

P. Huber, M. Lindner, and W. Winter, Comput. Phys. Commun. 167, 195 (2005), 
|hep-ph/0407333[ 

[81] M. Rolinec, Acta Phys. Polon. B37, 2049 (2006). 



[82] GLoBES, manual (2004), http://www.ph.tum.de/$~\sim$globes 



22 



[83] G. L. Fogli, E. Lisi, A. Marrone, and D. Montanino, Phys. Rev. D67, 093006 (2003), 
|hep-ph/0303064[ 

[84] M. Maltoni, T. Schwetz, M. A. Tortola, and J. W. F. Valle, New J. Phys. 6, 122 (2004), 



hep- ph/0405172| 



[85] J. N. Bahcall, M. C. Gonzalez- Garcia, and C. Pena-Garay, JHEP 08, 016 (2004), 
|hep-ph/0406294[ 

[86] A. Bandyopadhyay, S. Choubey, S. Goswami, S. T. Petcov, and D. P. Roy (2004), 
|hep-ph/0406328[ 

G. L. Fogli and E. Lisi, Phys. Rev. D54, 3667 (1996), [hep 1 ph/9604415| 

J. Burguet-Castell, M. B. Gavela, J. J. Gomez-Cadenas, P. Hernandez, and O. Mena, 
Nucl. Phys. B608, 301 (2001), [h^p" r ph/0103258[ 

[89] A. Donini, D. Meloni, and P. Migliozzi, Nucl. Phys. B646, 321 (2002), 



hep-ph/0206034 



[90] D. Autiero et a/., Eur. Phys. J. C33, 243 (2004), hep-ph/0305185 



[91] P. Lipari, Phys. Rev. D61, 113004 (2000), [hep r ph/9903481 , 



[92] P. Huber and W. Winter, Phys. Rev. D68, 037301 (2003), hep-ph/0301257 



23 



